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ABSTRACT 



During June and July, 198S, an AX13T survey was conducted southeast of Mindanao, 
Philippines in conjunction with CTD sections completed along the coast of Mindanao. 
Simultaneously, six Lagrangian drifters are launched in the Mindanao Current. Analysis 
of the data indicate a highly dynamic and complicated circulation which forms the origin 
of the North Equatorial Countercurrent. Strong shears arc observed across the 
Mindanao current and across a current located on the northeast coast of Borneo, and 
based on water mass characteristics and Lagrangian data, the general positions and 
strengths of the Mindanao and llalmahera Eddies. 



-y'(s> 3^ 

C.l 

TABLE OF COiNTLM S 

I. INTRODUCTION 1 

II. BATHYMETRY 5 

III. LAGRANGIAN DRIFTER MEASUREMENTS 8 

A. MINI-TRISTAR DRIFTER 8 

B. DRIFTER MEASUREMEN TS .- 8 

C. LAGRANGIAN TIUXJECTORIES 9 

1. Buoys launched along 7°N 11 

a. Buoy 52 II 

b. Buoy 51 13 

c. Buoy 55 13 

2. Buoys launched along 8°N 18 

a. Buoy 50 18 

b. Buoy 54 20 

c. Buoy 53 20 

D. LAGIUANGIAN VELOCITY CORRELATIONS 21 

1. Philippine Sea 23 

2. Celebes Sea 25 

E. LAGR.\NGIAN TEMPER/\TURES 26 

F. SUMMARY 30 

IV. AXBT AND CTD MEASUREMENTS 38 

A. VERTICAL PROFILES 38 

B. VERTICAL SECTIONS 39 

C. HORIZONTAL TEMPEIUM URE S I RUCTURE 42 

D. DYNAMIC STRUCTURE 43 

V. SUMMARY AND RECO.M.MENDATIONS 50 

A. SUM. MARY 50 

B. RECOMMENDATIONS 51 

iv 



‘’on?» 

APPENDIX A. AXBT TEMPEFLATLRE PROFILES 52 

APPENDIX B. LAG10\NGIAN TEMP1;RAH Ri; TIMF. Sl Rll.S 110 

LIST OF REFERENCES 124 

INITIAL DISTRIBUTION LIST 125 



LIST OF TAULtS 



Table 1. DRIFTER XLCROLOGV 13 

Table 2. TIME AVERAGED LAGRAXGIAN MEAN VELOCITIES 37 

Table 3. LAGRAXGIAN VARIANCES AND EDDY KINETIC ENERGY ... 37 
Table 4. DYNAMIC HEIGHT COMPARISONS 43 



LIST OF FIG LI RES 



Figure 1. Bathymeiiic Chart of the Survey Region 7 

Figure 2. MINI-TRISTAR Drogue and Surface Float Configuration 10 

Figure 3. Schematic of the Observed Motion of TRISTAR in Surface Waves ... 11 

Figure 4. Geometry’ of Location Calculations 12 

Figure 5. Lagrangian Trajectory of Buoy 52 14 

Figure 6. Buoy 52 Velocity Time Series 15 

Figure 7. Lagrangian Trajectory of Buoy 51 16 

Figure S. Buoy 51 Velocity Time Scries 17 

Figure 9. Lagrangian Trajectory of Buoy 55 IS 

Figure 10. Buoy 55 Velocity Time Scries 19 

Figure 11. Lagrangian Trajectory of Buoy 50 21 

Figure 12. Buoy 50 Velocity Time Series 22 

Figure 13. Lagrangian Trajectory of Buoy 54 23 

Figure 14. Buoy 54 Velocity Time Series 24 

Figure 15. Lagrangian Trajectory of Buoy 53 25 

Figure 16. Buoy 53 Velocity Time Scries Day 206-270 26 

Figure 17. Lagrangian Trajectoi-y of Buoy 53 27 

Figure IS. Buoy 53 Velocity Time Series Day 270-350 2S 

Figure 19.. Buoy 53 Velocity Time Series Day 340-430 29 

Figure 20. Buoy 51 Lagrangian N'elocity Correlations 30 

Figure 21. Buoy 52 Lagrangian Velocity Correlations 31 

Figure 22. Buoy 52 Lagrangian Velocity Correlations 32 

Figure 23. Buoy 53 Lagrangian Velocity Correlations 33 

Figure 24. Buoy 53 Lagrangian Velocity Correlations 33 

Figure 25. Buoy 50 La.srangian Velocity Correlations 34 

Figure 26. Buoy 54 Lagrangian Velocity Correlations 34 

Figure 27. Buoy 55 Lagrangian Velocity Correlations 35 

Figure 2S. Lagrangian Daily Minimum Sea Surface Temperature 36 

Figure 29. AXBT and CTD Station Locations 39 

Figure 30. Vertical Profiles of ‘femperature, Salinity, and Density .‘\nomaly 40 

Figure 31. Vertical Profile of Temperature, Salinity, and Density Anomaly 41 

vii 



Figure 32. Vertical Temperature Cross-section A 44 

Figure 33. Vertical Temperature Cross-section 13 45 

Figure 34. Vertical Temperature Cross-section (' 4o 

Figure 35. 23° Isotlicrm 47 

Figure 36. 14° Isotherm 47 

Figure 37. 14°-23° Thickness 4S 

Figure 3S. Buoy Trajectory composite 48 

Figure 39. Dynamic Heights 0, 300 m 49 



viii 



1. INIKODUCTION 



On Januar\‘ I, 19S5, ihe international TOGA (Tropical Ocean Global Atmosphere) 
program comtnenced. Planning for this program had been on going since 19S1 with its 
central focus being air-sea interaction and its elfect on the global climate. In 1982-1983, 
a strong El Nino event stimulated higher levels of interest in this research. TOGA is 
envisioned to be a 10-year program with its efforts centered around four major program 
elements. The first of these elements is modeling with an aim at fixing the predictability 
of the atmosphere through coupled air-ocean models. Second, empirical studies, based 
on historical data sets, will continue to contribute to the understanding of the processes 
that control the predictability of the environment. Third, process studies will be focused 
on those physical mechanism which are identified to be of critical importance to the 
understanding of large scale atmospheric - oceanographic variability. Finally, long term 
observations are necessary to determine the time-dependent structure of the tropical 
oceans and the global atmosphere. 

Although, global elTects of tropical air-sea interactions is central to the understand- 
ing of the circulation of the oceans at low latitudes along the western boundaries, yet 
until recently there had been few modern oceanographic observations in the western 
equatorial Pacific. Then, between 1985 and 1988, under the umbrella of the process 
studies element of TOGA, several expeditions were carried out as part of a Western 
Equatorial Pacific Ocean Circulation Study (WEPOCS). The third of these expeditions 
was a survey conducted by the R V Mocina If'aic called WEPOCS 111. The purpose of 
this cruise was [Lukas, et al, 1987]: 

1. defining the circulation of the low-latitude western boundaiy currents of the west- 
ern Pacific, including the possible cross-equatorial transports, 

2. determining the convergence of waters of northern and southern hemisphere origin 
in the far western equatorial Pacific, 

3. determining the source water mass c’ aracteristics for the Pacific to Indian Ocean 
Indonesian throughfiow, 

4. defining the Origin of the Equatorial L'ndercurrent, and 

5. determining the circulation patterns which supply the North Equatorial Counter- 
current near its origin. 

In order to learn more about the mesoscale structure in the region, the Naval 
Postgraduate School (NPS) supplemented the WEPOCS 111 program by providing ad- 
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ditional drifting buoys and an aircraft survey. Mini-TRISTAR drifting buoys were de- 
ployed to bracket the Mindanao Current. During three days in July, 1988, an Airborne 
Expendable Bathythermograph (.WBl') .survey was conducted by NFS southeast of 
Mindanao, in which 84 temperature proliles were obtained. 1 he data from these surveys 
are used in this study. 

Prior to TOGA and WEPOCS, our understanding of the circulation of the region 
had been based on ship drift reports and inJ'erences based on water mass properties 
[Wyrtki, 1961 and Wyrtki et ah, 1976]. The precision and accuracy of ship drift estimates 
is well known to be limited. Errors in ship positioning also limit the accuracy of these 
estimates. Wyrtki notes (1961) "The individual observation is an average over a period 
of 1 day and over a distance of about 400 km, somewhat larger than the typical 
midocean eddy, or the typical width of strong currents. The current speeds reported are 
consequently biased toward lower values. Included in the observed drift is also an un- 
determined efi'ect of wind on the ship, which largely depends on the loading character- 
istics of the particular ship," 

There have been several research cruises to the region. The sections made by the 
Kagoshi-maru (1952) and the Keiien-maru (1954-1956) were very effective at resolving a 
cyclonic vortex east of .Mindanao, called the .Mindanao Eddy [Takahashi, 1959]. Sections 
made by the Ailaniis II, at approximately 7°N and 8°N, and the survey conducted by 
the Takuyo (1965-1966) were effective in resolving the structure of the .Mindanao Cur- 
rent. From an analysis of this data. Cannon [Cannon, G. A., 1965] concludes "The off- 
shore limit of continuous strong currents is 70-80 km. Geosirophic surface currents and 
the ships set are consistent near Mindanao. The transport and the depth of the currents 
are diflerent on the two sections, which were only about 100 km apart." Since these 
expeditions took place before mesoscale eddies had been "discovered", they were not 
designed to adequately resolve eddies. These earlier investigations did not make direct 
measurements of ocean currents. 

The general circulation of the region has been described by Wyrtki [Wyrtki, 1961], 
The description that follows is for Jui and is based on his monograph. The North 
Equatorial Current (NEC) approaches ihe northeast coast of Mindanao from the cast 
and bifuricates into a northern flowing branch and a southern flowing branch. The 
southern branch, running from approximately 10°N to 3°N, is the Mindanao Current. 
The Mindanao Current flows along the coast of .Mindanao, and at its southern tip, also 
splits into two branches. One branch Hows southwestward over the ridge connecting 
Celebes and Mindanao (Celebes-Mindanao Ridge) while the second branch continues 



southward along the ridge. The South Equatorial Current (SEC) also approaches from 
the east along the northern New Guinea coast and at the northwest tip of New Guinea 
separates into two branches. The northern hianch then retrollecis into the North 
Equatorial Countercurrent (NECCj. 'Ihe second branch continues on between New 
Guinea and llalmahera until it reaches the Island ol' Ceram where a second bifurcation 
takes place, with one branch entering the Banda Sea and the other llowing northward 
into the Molucca Sea. This latter branch eventually enters the Pacific Ocean through 
the Molucca Passage, converges with the southward flowing Mindanao Current, and 
retroflects into the NECC. 

The southwest flowing extension of the Mindanao current which enters the Celebes 
Sea diverges into three separate extensions. The northern extension is very broad and 
passes through the Sulu Archipelago and into the Sulu Sea. The main extension passes 
through the central Celebes Basin and into the Java Sea via the Macassar Strait. The 
third extension is a retrollection of the current along the Celebes Peninsula which passes 
back over the ridge connecting Celebes and Mindanao between Celebes and the Island 
of Sangihe. This extension Hows into the convergence zone of the Mindanao Current 
and the SEC. 

Using water mass characteristics, and circulation patterns, Wyrtki [Wyrtki, 1961] 
suggested the presence of semi-permanent eddy features east of Mindanao (.Vlindanao 
eddy) and south of Mindanao (llalmahera Eddy). The exact size of these features is not 
well resolved although the radius is expected to be consistent with the Rossby radius 
measured in the region of 400 km [Emery, et al., 19S4J. The Mindanao eddy is located 
east of Mindanao at approximately 7°N, 130°E and is associated with the reversal of 
flow of the North Equatorial Current (NEC), the Mindanao Current and the NECC. 
The llalmahera Eddy is located at 3°N, 133°E and is associated with the convergence 
of the SEC and the NECC. 

This is a study of the physical oceanography of the region southeast of Mindanao 
and the flow of the Mindanao current into the Celebes Sea. The strategic importance 
of this region has been well documented and the effects of • lesoscale eddies on anti- 
submarine warfare cannot be overstated, for these reasons the circulation character- 
istics described in this study will be important to naval operations. The basis for this 
study will be the six TRISTAR Lagrangian buoys launched in the \icinity of the 
Mindanao Current, and the AXB f and CTD data collected in the region by the R/V 
Moana iVave. This study will be primarily descriptive in nature and limited to the sum- 
mer regime, as revealed by these data. Horizontal and vertical temperature character- 



3 



istics for the June-July (southwest monsoon) time frame for the region southeast of 
Mindanao will be discussed. Also included will be circulation patterns as determined by 
drifting buoys for the entire region from east of Mindanao to tlie Macassar Strait. 
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II. BAIIlYMCTin 



This study involves the throughllow of water between basins and would not be 
complete without a review of the bathymetrv’ that helps to govern the flow. Although 
Wyrtki [Wyrtki, 1961] provides a very thorough disscussion on the bathymetry of the 
region of this study, his discussion includes all southeast Asian waters. Since 
bathymetry', topography play such an important role in the magnitude and direction of 
ocean currents and eddies, a more detailed discussion is provided here. To avoid nu- 
merous references to Wyrtki, quotations will be used. Fig. 1, adapted from Wynki 
(1961) is a depiction of the bathymetry of the region. 

South of Mindanao is the Molucca Basin. Its boundaries are defined by a series of 
ridges running between Mindanao, 1 lalniahera. New Guinea, and Celebes. "The bottom 
of the basin is divided into five smaller basins which have depths between 3400 ni and 
4S00 m." The Molucca Basin is closed on the east by a ridge running between 
llalmahera and Mindanao with a sill depth of 2340 m. To the south, the boundary of 
the basin is defined by a ridge running between Celebes and New Guinea containing the 
islands of Taliabu, Mangolc and Obi. "The deepest threshold through this ridge is the 
Lifamatola Strait (located between Mangole and Obi) with a sill depth of ISOO m, which 
governs the replacement of the bottom water of all basins south of it." Completing the 
closure of the Molucca Basin is a ridge connecting the Celebes Peninsula with the 
southern tip of Mindanao. This ridge includes the islands of Siau and Sangihe and has 
a controlling sill depth of around 1400 m. 

Further to the west of the Molucca Basin is the Celebes Basin which has as its 
eastern boundary the ridge connecting Mindanao and Celebes. This basin is relatively 
flat and has a ma.ximum depth of 6220 m. The northern boundary of the Celebes basin 
is the Sulu Archipelago; "this is very shallow in large areas and consists of numerous 
coral banks." "Two channels, one east of Jolo with a sill of about 200 m depth, a-'d an- 
other east of Sibutu, with a sill of about 270 m depth, form slightly deeper but very- 
narrow thresholds." The western bouiidary of the Celebes Basin is the east coast of 
Borneo and the closure is completed in the south by the Celebes Peninsula. 

Providing for the passage of water out of the Celebes Sea, to the south, is the 
.Macassar Strait. Located between Borneo and Celebes, the central trench through this 
passage is around 2300 m in depth. 



The Philippine Basin, located to the east of Mindanao and to the north of New 
Guinea, is the region of the Mindanao eddy. Separated from the New Caroline Basin to 
the south by a ridge runnitw from IKdmahera through tltc Palau Islands to 't ap, its 
deepest depths of 10,000 m are found in the Philippine trench. Located at approximately 
S-S^N, 131°E is a seamount which shoals to a depth of 3000 m. 



fa 




Figure 1. Bathymetric Chart of the Survey Region: (adapted from Wyrtki, 1961) 
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111. LAGRANGIAN DRIFTER MEASUREMEN IS 



A. MINI-TRISTAR DRIFTER 

The data used in this study is a combination of temporal and spatial data and, as 
mentioned earlier, consists of drillers, AXBTs and CTDs. Here wo discuss the measure- 
ments made by the surface drilling buoys. 

Since a measurement of surface currents was desired, it was important to use a sur- 
face drifting buoy (subsequently called a "drifter") which would closely couple to the 
surface waters and minimize the effects of wind and waves. The buoy chosen is called a 
Mini-Tristar (Fig. 2) and was manufactured by Techocean Associates. The drifters fea- 
ture small spherical surface and subsurface floats and a large syinctrical drogue in the 
shape of a corner-radar reflector [Poulain et al.,iyS7]. The surface lloat was designed to 
remain on the surface at least SO'l'o of the time, while the subsurface float was designed 
with only weak negative buoyancy. This combination allowed slack to develop in the 
cable which connected these two floats. This would allow the drogue to move horizon- 
tally through the water with the vertical motion of the surface float decoupled from the 
horizontal subsurface motion. However, under strong wave forcing the surface float was 
observed to submerge (both in field tests and in a laboratory wave tank). The resulting 
buoyant force of the surface float would "take up" the slack in the line causing a dra- 
matic vertical oscillation of the drogue (Fig. 3), as the floats become coupled. Fortu- 
nately, the slippage is minimal. A detailed discussion of the characteristics and design 
of the TRISTAR drifters is given by Niiler (1987) [Niiler, et al., 1987]. 

B. DRIFTER MEASUREMENTS 

The buoys were tracked from two polar-orbiting NOAA satellites, using a system 
called ARGOS. The ARGOS system has been described by Poulain. [Poulain et al., 
1987, p. 3]. "Positions were determined from Doppler-shifted radio signals. The field of 
possible positions for a buoy is in the form of a half-cone (Fig.4) with the satellite at its 
apex, the satellite velocity vector as a.xis of symmetry, and the apex half-angle (A) such 
that: 

cos(A) = [(fr-fe);fej x c/V 



where: 

c = speed of light 



V = satellite speed relative to buoy 
fe = transmission frequency, and 
fr = receive frequency. 

The intersections of the various location cones with the ocean surface gives two possible 
positions of the buoy. Additional information - previous locations, range of possible 
speeds - is required to find which of the two positions is correct." 

The location quality of the fixes received from service ARGOS was based on the 
length of the satellite pass, the control over oscillator stability, geometric conditions, and 
how rapid the convergence of the fix was in the least-squares computation. The condi- 
tions encountered during tliis study provided, in most cases, a 6S% confidence interval 
that the fix was accurate to within 350 m. In many instances the fixes were accurate to 
within 150 m. Accordingly, the position data required only manual smoothing by re- 
moving obviously bad fixes. This was done by plotting the trajectories using the raw data 
and then eliminating those fixes which were clearly in error. In one instance. Buoy 53, 
a simple 5 point triangular filter was used to adequately resolve the initial recirculation 
it underwent in the .Mindanao liddy during the first month of the record. Additionally, 
due to problems with the buoy transmit terininal oscillator, approximately eight days 
of record were lost in the early part of Buoy 51's record. The temperature and 
submergence time of the buoys was also collected and were filtered using a cosine run- 
ning mean filter. 

Initial processing of drifter data was done by ser\ice ARGOS. Ihey provided 
monthly data tapes which contained location data as well as temperature and 
submergence records. The buoys in the study averaged 5.1 fixes per day. For each of the 
location fixes the quality was determined by service ARGOS. To obtain positions at 
standard times a linear interpolation was used on the position data to obtain positions 
ever>' twelve hours. Velocity time series were obtained from these interpolated positions 
and then smoothed using a 2.5 day running mean filter. 

C. LAGRANGIAN TRAJECTORIES 

The buoys were launched on 16 July and between 22 and 23 July, 19S8. The longest 
and shortest records were 209 and 17 days, respectively, with an average record length 
of 70 days. In order to obtain measurements of the shear across the .Mindanao Current, 
buoys were launched along lines perpendicular to the current at S°N and 7®N, with the 
position of the current determined by shipboard measurements. Table 1 lists the exact 
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Ministar Laaranoian Drifter 




Figure 2. MINI-TRISTAR Drogue and Surface Float Configuration 
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Figure 3. Schematic of the Observed Motion of TRISTAR in Surface 
Waves: (from Niiler, et al., 1987) 

dates and locations of launch in addition to the record length of each buoy and their 
failure mode. 

1. Buoys launched along 7°N 
a. Buoy 52 

Buoy 52 Mas launched on 16 July about 17 kni east of Tugubun Pt. on the 
south east corner of Mindanao in approximately 45U(J m of water (Figs. 5 and 6). Its 
initial irajcctor\' M'as to the southwest at a speed of 100 era s. By 19 July (Julian 200) it 
had slowed to 40 cm's as it began to pass over the Celebes-Mindanao Ridge. In moving 
over the ridge, the water depth shallows to a minimum of ISO m. Then, on 20 July, it 
backs to the southeast, nearly running aground on the small island of P Kawalusu lo- 
cated at 4.3°N, 125.4°E. On 22 July the buoy v\as established in a recirculation pattern 
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S,2-SUB satellite PO 
POSSIBLE POSITIONS 



Figure -1. Geometry of Location Calculations: (from Poulain et al., 1987) 



to the northeast and passed 22 km northeast of Sangihe througli 150 m of water. The 
buoy then continued toward the island of P Karakelong, then rapidly backed further to 
the north until it cleared the island on 28 July (Julian 209). Upon island passage, the 
buoy veered eastward and continued its drift into the NECC, and then accelerated to 120 
cm;s, the largest speeds observed by this which were the largest speeds observed by this 
buoy. It then recurved to the southeast, slowing to 115 cm s. On 1 August (Julian 213) 
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Table 1. DRIFTER NECROLOGY 



BUOY 


Date of 
launch 


Location of 
launch 


Record 

length 


1 ailure 
mode 


Buoy 50 


23 July 88 


7.958°N 

126.556“E 


48 days 


Lost 

Drogue 


Buoy 51 


16 July 88 


6.983“N 

127.304“E 


33 days 


Lost 

Drogue 


Buoy 52 


16 July 88 


6.944“X 

126.619“E 


53 days 


Lost 

Drogue 


Buoy 53 


22 July 88 


7.983°N 

127.95SM: 


220 days 


Grounded 


Buoy 54 


22 July 88 


7.990“N 

126.975°E 


17 days 


Picked up 
at Sea 


Buoy 55 


16 July 88 


6.86rX 

126.86rE 


63 days 


Grounded 



it entered an anticydonic oscillation in the Philippine Sea. When the buoy reached 
7°N, 13-4°E it slowed to 60 cm s and began a cyclonic recirculatio.:. It completed a 75 
km radius gyre, maintaining 60 cm's, by 23 August (Julian 235) and then began a second 
anticydonic trajectory. A quarter of a period into this oscillation the buoy is lost on 17 
September (Julian 260). 

b. Buoy 5 1 

■ Buoy 51 (Figs. 7 and 8) was launched on 16 July about 75 km east of Buoy 
50. Its initial trajectory was to the south at 44 cm, s. On 18 July it began to slow to 20 
cm, s. At this point, the buoy transmit terminal oscillator I'ailed. I he oscillator began 
functioning again on 2S July (Julian 209), and the buoy was located at 6.7°N, 129.4°E. 
The buoy began circulating to the south and east in a cyclonic pattern at 66 cm s. It 
continued in this pattern until 19 August (Julian 231) when it slowed to 7 cm, s before 
sharply ti -ning back to the west. The buoy was last tracked at 6.3°N, 133.6°E. It was 
lost on 20 August 8S (Julian date 232). 

c. Buoy 55 

Buoy 55 was launched on the axis of the Mindanao Current on 16 July and 
iinmediately began to drift to the southwest at a speed of 90 cm s. Within 24 hrs it had 
passed within 5000 m of Cape Augustine and began to pass over the Celebes-Mindanao 
Ridge wherein it bridly maintained a speed of 80 cm s (Figs. 9 and 10). By 22 July it 
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Figure 5. Lagrangian Trajectory of Buoy 52: Julian dates mark the beginning and 

end of each record and a square is placed as a marker for 12z every sec- 
ond day. 

had cleared the ridge and began to slow. During one 24- hour period between 24 and 25 
July, the speed dropped to as low as 2 cm's. On 26 July the speed increased to 25 cm's 
and the buoy began to drift back to the south, eventually completely retrollected back 
to the east and then drifted along the coast of the Celebes Peninsula, passi:., within 5000 
m of landfall. It began to move back to the north on 15 August (Julian 227) and, after 
two days, along this track it turned sharply to the west. During this evolution, the buoy- 
maintained a speed of 25-40 cm s with the exception of a brief period on 17 August 
(Julian 229) when it slowed to 12 cm's. By 19 August the buoy had completed an oscil- 
lation with a radius of 150 km, had accelerated to 80 cm/s and was moving toward the 
southwest extreme of the Sulu Archipelago. On 23 August (Julian 235) it passed within 
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Figure 6. Buoy 52 \'elocity Time Series: N'elocities were calculated from in- 

terpolated trajectories and then filtered using a 2.5 day triangular run- 
ning mean. The dashed line is the v component, the solid line is the u 
component and the chain-dot line is the modulu.s. 

10 km of the Subutu Passage through the Archipelago but managed to maintain its 
trajectory around the Celebes basin. It eventually moved back to the south within 10 
km of the coast of Borneo. On 31 August the buoy entered the .Macassar Strait. On the 
day prior to entering the Strait, the buoy had reached its maximum sustained speed, 1-45 
cm/s. .Moving through the Macassar Strait, the buoy meandered along the 1400 m 
isobath with a wavelength of 500 km. During the last 5 days in which the buoy transited 
the Strait, it maintained a consistent 35 cm s speed until it apparently lost its drogue. 
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Figure 7. Lagrangian Trajectory of Buoy 51: Julian dates mark the beginning and 
end of each record and a square is placed as a marker for I2z ever>’ sec- 
ond day. 

having either run aground or been trapped in fishing nets, on 18 September, as it was 
entering the Flores Sea. 
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Figure 8. Buoy 51 \’elocity Time Series: \’elocitics were calculated from in- 

terpolated trajectories and then filtered using a 2.5 day triangular run- 
ning mean. The dashed line is the v component, the solid line is the u 
component and the chain-dot line is the modulus. Gap in record rellects 
failure of buoy transmit terminal oscillator. 
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Figure 9. Lagrangian Trajectory of Buoy 55: Julian dates mark the beginning and 
end of each record and a square is placed as a marker for 12z ever} sec- 
ond day. 



2. Buoys launched along 8°N 
a. Buoy 50 

Buoy 50 was launched on 23 July, 14 km northeast of Bingai Pt. on 
Mindanao and moved to the sout rt 120 cm's on the inshore side of the Mindanao 
Current. During its first day of drifi, the buoy passed within 4000 m of Pusan Point and 
then continued to the south (Pigs. 11 and 12). On 24 July, at appro.ximately 6.9°N, 
126.5°E, the buoy veered to the southeast and slowed slightly. On 25 July, the buoy be- 
gan to drift farther to the southeast and passed Cape San Augustine 3500 m to the 
southwest. Continuing to the southwest, the speed varied from 40 to 50 cm/s as it 
passed over the Celebes-.Mindanao Ridge at Sarangani Island in ISO m of water. Por 
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Figure 10. Buoy 55 \'elocity Time Series: Velocities \\ere calculated from in- 

terpolated trajectories and then filtered using a 2.5 day triangular run- 
ning mean. The dashed line is the v component, the solid line is the u 
component and the chain-dot line is the modulus. 

the next few days the buoy continued to move to the southwest with little variability in 
speed. Then on 3 August, the buoy slowed to 10 cm s as it entered a small cyclonic os- 
cillation at 3°N, 122. 5°E. Upon completing this rotation, the buoy moved to the north- 
west at 30 cm s on 8 August. The buoy approached wihin 100 niilcs of the Sulu 
Archipelago, and then turned to the southeast, following the 4000 m isobath. 

As the buoy approached the entrance to the Macassar Strait, its velocity 
accelerated to 80 cm s. It passed the entrance on 30 August (Julian 242) and then passed 
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within 2000 m of the northern tip of Celebes wherein it turned northward into the Cen- 
tral Celebes Sea. Near 2°N, 121 °h the buoy speed was 5 cm s. d'he speed brielly in- 
creased to 16 cm's two days prior to the buoy being Iom on 1 1 Sepieinber (.lulian 2>A). 

b. Buoy 54 

Buoy 54 was launched on 22 July on the axis of the .Mindanao Current, as 
determined by shipboard measurements (Figs. 13 and 14). It initially moved to the south 
with a speed of 80 cm;s. An interesting aspect of this buoy’s trajectory can be made by 
comparing its movement with the movement of Buoy 52, which was launched 7 days 
earlier and to the southeast of Buoy 54's launch site. The relative positions of It was 
expected that , based on the relative positions of their two tracks. Buoy 54 to the north 
and west of Buoy 52, Buoy 54 would retroflect back to the northeast in a manner similar 
to Buoy 52. Contrarily, Buoy 54 continued over the Celebes-Mindanao Ridge at a speed 
between 50 and 70 cm,'s. In the process of passing over the ridge, the water depth shoals 
to 150 m and the buoy passed within 1000 m of the small islands of P Kawio and P 
Kamboling, located 90 km north of Sangihe. By 25 July, the buoy accelerated to 100 
cm;s and is continued along a path similar to the path of Buoy 50. When the buoy 
reached the central Celebes Sea, it had slowed to 50 cm's, it was beginning to move to 
the south when it w’as picked up on 10 August (Julian 222) by a ship and taken to the 
southern coast of .Mindanao. 

c. Buoy 53 

Buoy 53 was launched on 22 July, east of the .Mindanao Current. This buoy 
has the longest record (220 days) and spent its entire time in the Philippine Sea (Figs. 
15 through 19). Initially, it moved in a slow elliptical path, with the major and minor 
axis of this circulation 100 km and 50 km, respectively. It remained in this circulation for 
the first 42 days of its history mox’ing at speeds between 10 and 40 cm s with occasional 
accelerations to speeds as high as 70 cm,s. The center of this circulation was 7°N, 
12S.5°E and most likely captures the signature of the center of the .Mindanao Eddy. On 
2 September (Julian 245) the buoy began to move to the southeast and slowly acceler- 
ated to 80 cm's . 22 Sept, at which time it moved eastward in the NECC. With its 
speed as high as 50 cm.'s to 100 cm;s, the buoy continued to the east reaching 141. 7°E 
on 8 October. At this point it moved to the north, and then northwesterly. On 7 No- 
vember (Julian 31 1) it had completed a retrollection to the west. It then traced out a 400 
km wavelength oscillation to the north, with typical speeds of 25 to 30 cm/s, passing 200 
km east of Palau, and entering the flow of the NEC. Once in the NEC it moved 
westerly, passing through a 650 km wavelength oscillation. Speeds in the NEC varied 
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Figure 11. Lagrangian Trajectory of Buoy 50: Julian dates mark the beginning 
and end of each record and a square is placed as a marker for I2z every 
second day. 

between 5 and 60 cm s with the largest speeds at ITN, 12S°E. The buoy apparently 
drifted ashore or lost its drogue on the Island of Luzon on 26 Lcbruary. 

D. LAGRANGIAN VELOCITY CORREL.\TIONS 

The velocity time series were linearly interpolated using a great circle solution. 
Similar to the method used by Poulain [Poulain, et al, 1^87, pp. 4), the smooth (2.5-day 
running mean filter) velocity time series were used to calculate the mean velocities f^r 
each drifter. From the mean velocity, a residual velocity data set was computed and 
used to calculate the correlation functions. "Lagrangian time lagged velocity corre- 
lations were then computed by averaging time lagged velocity products from individual 
buoys." Using auto-correlation function and cross-correlation functions; 

R(iX = imr)iXi + r)j 
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Figure 12. Buoy 50 Velocity Time Series: Velocities were calculated from in- 

terpolated trajectories and then filtered using a 2.5 day triangular run- 
ning mean. The dashed line is the v component, the solid line is the u 
component and the chain-dot line is the modulus. 

^(t/l% = £lT-(/)r(; + T)J 
^(KF-)_, = £lLX/)r(r-T)] 

where t is the time lag. These functions are shown in Figs. 20 through 27 with redun- 
dant halves of the curves, due to symmetry, omitted. Using the velocity time series and 




Figure 13. Lagrangian Trajectory of Buoy 54: Julian dates mark the beginning 
and end of each record and a square is placed as a marker for 12z every 
second day. 

the correlation functions, Tables 2 and 3 were constructed to compare the mean veloci- 
ties, variances and eddy kinetic energy (EKE) per unit mass of each buoy. 

1. Philippine Sea 

Buoys 51, 52, and 53 remained in the Phillipine Sea. Their Lagrangian velocity 
correlations (Figs. 20 thru 24) show that most of the din'erences discernible between 
these correlations are dependent upon their launch position and initial movement. Buoy 
52 has the highest variance and energy associated with it, however; it was initially in the 
strong flow of the Mindanao Current. If the effects of the current are removed (Fig. 22), 
we see a correlation more closely resembling the other two buoys. (Figs. 2U and 21). 
This also leads to a velocity modulus of 4. S3 cm s which compare more closely to Buoys 
51 and 53. For all buoys in the Philippine Sea, the energy is concentrated in the zonal 



23 



160.0- 



120.0- 




time (DAYS) 



Figure 14. Buoy 54 Velocity Time Series: Velocities were calculated from in- 

terpolated trajectories and then filtered using a 2.5 day triangular run- 
ning mean. The dashed line is the v component, the solid line is the u 
component and the chain-dot line is the modulus. 

direction, even for Buoy 52 which was in the southwesterly flow of the Mindanao Cur- 
rent. 

The most significant difference in time scales is that of Buoy 53 , Fig. 23. The 
early part of this buoy's trajector>’ was trapped in a slow oscillation (probably the center 
of the .Mindanao Eddy), with a time scale 30 days. However, if the early part of the re- 
cord from this buoy's velocity time series is removed and the correlations (Fig. 24) re- 
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Figure 15. Lagrangian Trajectory of Buoy 53: Julian daies mark the beginning 
and end of each record and a square is placed as a marker for 12z every’ 
second day. 

computed, the result is a lime scale of the order of 10 days, uhich is exactly what was 
observed for Buoys 51 and 52. 

2. Celebes Sea 

Buoys 50, 5-1, and 55 all entered and remained in the Celebes Sea (Figs. 25 
through 27). The slightly lower kinetic energy of Buoy 50 is associated w’ith the small 
oscillation it undergoes between 3 August and 8 August. Regardless of this oscillation 
the total kinetic energy of the buoy is of the order of magnitude of Buoys 54 and 55. The 
energy of the buoys appears to be evenly distributed between the north-south and east- 
west planes with the exception of Buoy 54 which was picked up at sea just as it began 
a southward excursion. For all three cases in the Celebes Sea, the velocity time scales 
of ten days is the same as that found in the California Current system by Poulain 
(Poulain, et al., 1987). 
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Figure 16. Buoy 53 Velocity Time Series Day 206-27U: Velocities were calculated 

from interpolated trajectories and then filtered using a 2.5 day triangu- 
lar running mean. The dashed line is tl^e v component, the solid line 
is the u component and the chain-dot line is the modulus. 

E. LAGRANGIAN TEMPERATURES 

Sea surface temperature was measured inside the surface float by a thermistor im- 
bedded in a stainless steel bolt which passes through the shell of the float beneath the 
waterline. The measured temperature exhibited diurnal changes ranging from 2S°C to 
33“C (see Appendix B). It is not certain whether this diurnal variability is real or not. 
For example, on a calm day, the surface float might absorb more of the sun's radiation 
than the surrounding water. Since there was no way of compensating for dilfercntial 
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Figure 17. Lagrangian Trajectory of Buoy 53: First 53 Jays while in the 

Mindanao EJJy. 

heating of the surface float and water, and my main interests were the larger scale, lower 
frequency patterns, an attempt to remove the diurnal variability was made. This was 
done by using only the mean minimum daily temperature for each buoy. On average this 
would be a record of the mean nightime sea surface temperature and the diurnal elTects 
would be filtered out. These temperature reeorJsas a time series are presented in Fig. 28. 
It should be kept in mind that the record is also a spatial series. 

In general, the mean temperature variations between buoys is small. The temper- 
ature spread is usually between 29°C and 29.5°C especially between 7 August (Julian 
219} and 29 August (Julian 241 ). However, most buoys experienced anomalous high and 
low temperatures. For example, Buoy 55 (launched on the inshore side of the Mindanao 
current) had an initial temperature of 28.1° C, considerably lower than the other buoys 
launched on the inshore side. This could be explained by deep water upwelling along the 
coast due to the prevailing southwest monsoon winds. The large drop in temperature of 
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VELOCITY (CM/S) 




Figure 18. Buoy 53 ^'elocity Time Series Day 270-35U: Velocities were calculated 

from interpolated trajectories and then filtered using a 2.5 day triangu- 
lar running mean. 

Buoy 52 (launched inshore and north of Buoy 55} as it drifted to the south, on 3 August, 
(Julian 215) also suggests the presence of colder water along the shore. Between August 
16 (Julian 22S) and August 28, Buoy 55 e.xperienccd unusually high temperatures of 
29.5°C and 29.6°C. During these episodes, Buoy 55 also experienced brief but small ac- 
celerations in its speed. This suggests that the increased temperatures occur in con- 
junction with a strengthening of the temperature gradient in the Central Celebes Sea and 
along the coast of the Celebes Peninsula. On September S (Julian 251), as Buoy 55 is 
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Figure 19. Buoy 53 Velocity Time Series Day 3-4U-I30: (Julian days greater than 

366 arc in reference to 19S9). Velocities were calculated from interpo- 
lated trajectories and then filtered u^ing a 2.5 day triangular lunning 
mean. 

beginning to enter the Flores Sea, it experienced a rapid temperature decrease from 
29.4°C to 2S.6°C. 

Between 27 July (Julian 20S) and 1 August, ,Buoy 50 e.xperienced a marked increase 
in temperature from 29.5“C to 3U.d“C. At the same time. Buoy 50 accelerated Irom 30 
cm;s to 50 em s and was located at 4.2»N, 124.5“i:. A speculative explanation for tins 
marked increase in temperature and the concurrent acceleration could be that this Id- 



29 



(CM/S)**: 



2500.0-t 



2000-0- 

1500.0- 

1000.0- 
C\1 

500.0- 
0.0- 
-500.0 
-1000.0 
-1500.0 

-50.0 -30.0 -10.0 10.0 30.0 oO.O 

TIME LAG (DAYS) 

Figure 20. Buoy 51 Lagrangian \’elotity Correlations 

cation marks the boundary’ between the water mass of a south equatorial Pacific intru- 
sion and the resident Celebes Sea water mass. 

Buoys 52 and 53 both entered the NhCC after separating from the influence of the 
Mindanao and llalmahera Eddies around 24 August (Julian 236) and iimnediately ex- 
perienced a warming from 29.5°C to 3U°C 

showing that the core of the NECC is likely marked by a meridional temperature max- 
imum. In both cases this warming is concurrent with accelerations in speed. In the case 
of Buoy 52, the speed increases from 3U em s to SO cm s by 29 August. Buoy 53 experi- 
enced minor warming episodes on 4, 12, and 17 September and in each case there was 
a concurrent acceleration between 5 and 10 cm s. 

F. SUMMARY 

It is difficult to generalize circulation patterns from only six dn ^ers. Nevertheless, 
I will tr>’ to summarize these results. In relation to the .Mindanao Current, three drifters 
were launched on the inshore side, two were launched on the seaward side, and one was 
launched on the axis. The buoys launched inshore or on the axis all moved initially to- 
wards the Celebes Sea at speeds in excess of lOO cm.'s and with kinetic energy about 1000 
{cmjsy. While the buoys launched to the seaward side also had a high energy relative 
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(CM/S)** 




10 niidocean values, the energy was considerably less, (0(500)), than the Mindanao 
Current. With the exceptions of Buoys 50 and 55, both launched inshore , the energy 
was predominately associated with the zonal How. In the cases of Buoys 50 and 55, the 
energy was evenly distributed between zonal and meridional flow. The time-averaged 
mean velocities also reflect this distribution of energy. At the surface the Mindanao 
current was strong and narrow. 

The trajectories of Buoys 50 and 55 suggest that the Celebes Sea was dominated by 
a large cyclonic gyre in the central basin and a strong southeasterly current along the 
coast of Borneo which bifurcates at the entrance to the Macassar Strait. Buoy 55 which 
reaches a maximum sustained speed of 145 cm s passed within 50 km of the coast of 
Borneo and then enters the Strait. Buoy 50 reached a maximum sustained speed of 80 
cm, s, approached the coast of Borneo 60 km further to the east, and then bypassed the 
entrance to the Strait and moved northward in the Celebes Sea. Surlace How into the 
.Macassar Strait seemed to be restricted to waters close to shore and northwest of the 
Strait. 

The trajectories east of Mindanao reveal three distinct circulations. The .Mindanao 
eddy center was located at 7°N, 12S.5'’E, as indicated by the intitial trajectory of Buoy 
53. Buoy 52 completes a anticyclonic oscillation "over the top" of the llalmahera Eddy 
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Figure 22. Buoy 52 Lagrangian Velocity Correlations: Days 22U-25U. computed 

alter removing ell'ccts of the Mindanao Current 



and then makes an additional cyclonic circulation lurther to the east revealing the com- 
plexity of the circulation associated with the origin of the NECC. The strongest flow 
was encountered at 5°X, 12S°E. fhis is a region of a front between SEC and NEC de- 
rived waters and will be discussed in the next section. 
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Figure 23. Buoy 53 Lagrangiaii ^'elocity Correlations: Days 204-424 




Figure 24. Buoy 53 Lagrangiaii ^'elocity Correlations: Computed alter removing 

the elTects of the Mindanao Eddy. Days 204-257. 
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Figure 25. Buoy 5U Lagraiigiaii \’eIocity Correlations 




Figure 26. Buoy 54 Lagrangiaii Velocity Correlations 
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Figure 27. Buoy 55 Lagrangiaii \’elocity Correlations 
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Figure 28. Lagrangian Daily Minimum Sea Surface Temperature 
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Table 2. TIME A\ TRACED LAGRANClAN MEAN \ ELOCITIES 



BUOY 


Mean L 


.Moan V 


.Mean 




velocity 


N'clociiy 


Speed 




(cm s) 


(cm s) 


(cm s) 


Buoy 50 


-13.9 


-9.S 


17.U 


Buoy 51 


5.4 


-7.4 


14.6 


Buoy 52 


29.4 


2.2 


29.5 


Buoy 53 


6.6 


3.0 


7.3 


Buoy 54 


-25.7 


-16.4 


30.5 


Buoy 55 


-18.6 


-22.0 


28.8 



Table 3. LaGRaNGIAN VARIANCES AND EDDY KINETIC ENERGY 



Buoy 


L Variance 


V variance 


U EKE 


V EKE 




{cmlsf 




{cmisf 


{anisf 


Buoy 50 


920 


830 


4o0 


415 


Buoy 51 


1590 


816 


795 


408 


Buoy 52 


1178 


1U02 


589 


501 


Buoy 53 


874 


196 


437 


98 


Buoy 54 


1269 


792 


635 


396 


Buoy 55 


1470 


1476 


735 


738 
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IV. AXI3T AND CTD MKASUKLMLN l b 



The AXBT survey conducted by the Naval Postgraduate School in July, 198S was 
confined to a region southeast of.Vlindanao. The zonal and meridional boundaries of the 
survey were PN-T^N and 126“E-1 30°l£, respectively (Tig. 29). Areal coverage was lim- 
ited by available flight time and equipment failure. The region was chosen to comple- 
ment the ship survey and to further document the mesoscale variability. The processing 
of the 85 profiles used in this study was conducted at the Naval Postgraduate School. 
While 95% of the data could be adequately smoothed by linearly interpolating across 
bad data points, the remaining 5% required a cosine running mean filter for smoothing. 

The Cl D casts used in this study were collected during the second leg of the 
WEPOCS III cruise. Casts were made to the bottom or 4500 m depending on the local 
water depth. The first section occupied by the R:V Moana ITave was from 7°N, 134®E 
to 2.5°N, 131.7°E, with stations spaced half a degree of latitude apart. The ship then 
proceeded to occupy stations along the southwest coast of Mindanao with the intention 
of mapping the coastal circulation. Stations were occupied parallel and perpendicular to 
the coastline from Moro Gulf to 7®N on the Pacific side of Mindanao (Fig. 29). .Most 
station spacing were 15 nm, near shore to Mindanao. 

Both the CTD sampling and processing are similar to that used during the WEPOCS 
II cruise and are described in the WEPOCS 11 data report. [Tsuchiya, et ah, 1987]. "The 
data were obtained by a modified Neil Brown Mark 111 CTD. The CTD data were ini- 
tially processed into a filtered 1 2 second average time-series during data acquisition. 
This time series was then pressure-sequenced into 2-deeibar pressure intervals. A spike 
filter was employed to remove large pressure, temperature and conductivity spikes from 
the time-series data and a "ship-roll" filter was applied to disallow pressure reversals." 
I will focus my discussion of these data to the thermal structure above 300 m. Vertical 
profiles, sections, and horizontal charts will be e.xamined to determine the geographic 
variability and flow patterns. 

A. VERTICAL PROFILES 

The general vertical structure of the region is represented by two CTD profiles (Figs. 
30 and 31) selected for comparison at opposite extremes of the survey area. One, near 
7°N, 127°E, station 70, acquired on 16 July, 1988, was in the vicinity of the Mindanao 
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Figure 29. AXBT and CTD Station Locations 

Current. The other, Station 32, located near 3°N, Ul^F, in the vicinity of the Halmahcra 
Eddy, was acquired on 8 July, 1988. 

The vertical temperature distributions at these stations show a notable din'erence in 
the depth and thickness of the therniocline. In the southern waters, the ihermocline is 
located between 100 m and 300 depth, whereas in the northern waters it is located be- 
tween 30 m and 150 m depth. The salinity maximum represents the Subtropical Lower 
Waters (Wyrtki, 1961]. The northern station has a salinity maximum of 35 ppt and it is 
located 50 m below the mixed layer at 100 m. The southern station has a slightly higher 
salinity maximum at 35.2 ppt and is located at 200 m or 100 m below the mixed layer. 

B. VERTICAL SECTIONS 

The vertical cross-section plotted at 6.5“N from 127® E to 12S.6®E {Fig. 32) , be- 
ginning 100 km from the southern coast of Mindanao, shows isotherms generally sloping 
down toward the coast, consistent with equatorial flow. The greatest slopes are below 
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Figure 30. Vertical Prol'iles of Temperature, Salinity, and Density Anomaly: CTD 

station 70 located at 6.95°X 126.7TE 

150 db between 127°E and 127.4°E, and strengthen with depth down to 3U0 m, which 
was the limit of the AXBT. Isotherms 22°C-2S°C, above 125 m, slope up toward the 
coast, evidence of possible upwelling. The 10°C-22°C isotherms slope down toward the 
east at 128. 3°E marking the location of the Mindanao eddy. Surface temperatures are 
about 30°C. Temperatures at 300 db are less than 10°C. 
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Figure 31. \’ertical Prolile of Temperature, Salinity, and Density Anomaly: CTD 

station 32 located at 2.S>r\ 131.39°T 

The north-south cross-section located along 129°F (Fig. 33), and the north-south 
cross-section located along 130®F (Fig. 34) are through the llalmahera Eddy. Both of 
these cross-sections have a weaker thcrmocline gradient than Fig, 32. However, toward 
the north there is a tightening of the gradient which more closely resembles that found 
off the coast of .Mindanao (Fig. 32). 'fhe ina.ximum surface temperature observed in the 
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two cross-sections is along 129°L-. and is 30°C. In both cross-sections there is a broad 
warm core feature between 2.5°N and 5.0°N, It extends to 300 db and the vertical dis- 
placement of istherms is greater wiili depth. In tlie western section ll ig. 3.1) there is a 
stronger gradient immediately below the mixed layer and a 300 db temperature of less 
than 12°C south of 3.3°X and greater than 12°C north of 3.3°N. The eastern section 
(Fig. 34) in addition to having a weaker gradient below the mixed layer, has considerably 
warmer deep water south of 3.0‘’N, 14°C and cooler deep water north of 3.8°N, 12°C. 
Surface temperatures greater than 30®C are formed only close to 5“N at 129°E. 

C. HORIZONTAL TEMPERATURE STRUCTURE 

I have chosen two isotherms to characterize the shape of the ihermocline in the re- 
gion, 23°C and 14°C.The 23° isotherm lies close to but below the mixed layer in these 
data. The 14°C isotherm depth is at the bottom of the thermociine. Its depth exceeded 
300 m in a few cases but was measured for all casts. 

The depth of the 23°C isotherm (Fig. 35) ranges from 76 m to 169 m. At S°N the 
depths increase to the coast (to about 140 m), marking the strong southwest flow of the 
Mindanao current. The region of shallow depths extends from the Mindanao eddy to the 
southwest and is outlined by the 110 m contour. In the llalmahera eddy to the south- 
east, in the area of the SEC, depths are about 130 m. The Mindanao eddy, located at 
approximately 1 28. 5°N, 7.0°E, is associated with the 23°C isotherm depths shallower 
than 90 m. To the south, across a frontal boundary which separates the NEC water and 
SEC water, the llalmahera eddy is located northeast of llalmahera at approximately 
129°N, 4°E. Associated with this eddy, the 23°C isotherm deepens to 160 m. 

The depth of the 14°C isotherm (Fig. 3o) ranges from 135 in to 320 m. Patterns are 
similar to those of the 23°C isotherm but are better defined and stronger. A 20 m con- 
tour inter\al is used here while a 10 ni contour interval is used with the 23°C isotherms. 
The 160 m contour at 7°N, 12S°E marks the location of the Mindanao eddy and is 
shifted slightly to the southwest of the location at 23°C. At 8°N, the 14°C isotherm 
deepens to 260 m, so the Mindanao current is strong and clearly delineated. Going to 
the southwest across the frontal boundary, the 14°C isotherm deepens from 160 m to 
240 m. The Halmahera eddy is marked by isotherm depths greater than 300 m. 

The depth difference between these twm isotherms is a measure of thermociine 
thickness and is shown in Fig. 37. The dilferences range from 50 to 20 m. 

The geographic distribution seems somewhat bivariate; waters from north of the 
equator have thicknesses of about 80 m while those from the southern hemisphere have 
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thicknesses exceeding 12U m. The boundary between these waters is marked by a strong 
gradient. The location of the Mindanao eddy is best indicated by the dosed CU m con- 
tour at 6°N, 128. 5°E which is slightly south oftlie location indicated on the other charts. 
The eddy olT llalmahera is marked by closed contours with thicknesses exceeding loU 
m. Fig. 38, a composite of all the buoy trajectories, is included here to show the motion 
of the buoys relative to the hydrographic features. The .Mindanao eddy is delineated east 
of Mindanao by the cyclonic circulation completed by Buoy 53. Although the 
Halmahera eddy does not appear on this figure, the anticyclonic recurving pattern across 
the Philippine Sea is consistent with the prescence of a warm core eddy north of 
Halmahera, The eddy located further to the east at 5“N, 134‘’E was not resolved with the 
AXBT data. 

D. DYNAMIC STRUCTURE 

In order to see if the dynamic structure in the region could be calculated from AXBT 
data, dynamic heights were calculated for the two CTD stations discussed above. First, 
dynamic height was calculated using the observed salinities and then they were again 
calculated using a constant salinity of 34.5 ppt. The results are shown in Table 4. From 
the dilTerence calculations in fable 4 a maximum error of 15To in the dynamic height 
calculation should be e.xpected when assuming constant salinity. This error was ac- 
ceptable for further calculations. 

Using the CTD data with observed salinities and the AXB f data and an assumed 
constant salinity of 34.5 ppt, a dynamic height field was calculated for the entire region 
(Fig. 39). The highest dynamic heights, 1.4 dyn m, are in the region to the south where 
the thickness is at a maximum. Also, the strong gradient can be seen in this field between 
the waters from the northern hemisphere and the southern hemisphere. The geostrophic 
velocity was calculated at this front to be about 70 em s, about half the velocity meas- 
ured by the drifters. In the northern region, dynamic heights are lower, 0.9 to 1.0 dyn 
m, and the Mindanao eddy does not have a strong signature. 



Table 4. DYNAMIC HEIGHT CO.MPARISONS 



Station 


Constant Salinity 
Dynamic Height 


Observed Salinity 
Dynamic Height 


32 


1.2815 


1.259 


70 


1.0568 


.98476 


difference 


.23 


.27 
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EAST LONGITUDES 



Figure 32. N'ertical Temperature Cross-section A: Located along 6.5°N 
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NORTH LATITUDES 

Figure 33. Vertical Temperature Cross-section B: Located along 129°H 
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Figure 3-4. \'ertical Temperature Cross-sectioii C: Located along IjO C 
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Figure 36. 14° Isotherm 
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Figure 38. Buoy Trajectory composite 
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Figure 39. D\namic Heights 0/300 m 
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V. SUMMARY AND RECOMMENDATIONS 



Prior to the commencement of TOGA and the subsequent WHPOCS cruises, our 
understanding of the circulations associated with the western equatorial Pacific and the 
origin of the NECC had been based upon ship drift reports and widely separated water 
mass samples. There had been no direct attempt at measuring the current speeds nor had 
a large scale effort been made to determine the dynamic structure based on closely 
spaced samples. Hopefully, this study has played a role in closing the gap between the 
heretofore inference and a true observed description of the circulations associated with 
the Mindanao and Halmahera Eddies. 

The major need for an increased understanding of the circulations in the western 
equatorial Pacific arises from its link to understanding the global climate (El Nino). The 
first link to prediction is the climatology of the region of concern. This study should 
serve as a beginning toward a better understanding of the sunmiertime circulation in the 
western equatorial Pacific. An understanding of these circulations also has tactical mil- 
itary significance to the United States. The largest forward deployed Naval bases are 
located in the Philippine Islands and the primary Battle Group Pacific- Indian Ocean 
transit route is through the western equatorial Pacific and the Indonesian Archipelago. 
Vloreover, every submarine olTicer knows the importance of being able to predict and 
locate mesoscale ocean eddies for Anti-Submarine Warefare. 

A. SUMMARY 

This study has been an investigation into the circulation associated with the 
Mindanao and Halmahera eddies in June and July, 19SS. The drifter data, as summa- 
rized earlier, showed a strong cyclonic shear across the Mindanao current as evidenced 
by the dificrence between the maximum velocity on the inshore side of 115 cm, s (Buoy 
50) and the minimum velocity on the seaward side of 20 cm s (Buoy 53). 'fhe buoys in 
the Philippine Sea were o' lived to interact with three dilferent eddies. The Mindanao 
Eddy located at 7°X, 12S.„°E, the Halmahera Eddy located at d°N, 129. 5°E, and a third 
unnamed eddy located at 5°N, 133°E. These Eddies mark the origin of the NECC and 
provide a mechanism for recirculation of water back into the NEC and SEC. 

The buoys that drifted into the Celebes Sea were all obser\ ed to drift around a large 
cyclonic gyre that dominated the central Celebes Sea. Buoys 50 and 55 both drifted close 
to the coast of Borneo and toward the entrance to the .Macassar Strait. However, at the 
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entrance to the Strait Buoy 50 retrollects buck to the north and Buoy 55 continues 
through the Strait. Both buoys experienced an increase in speed along the northeast 
coast of Borneo. Buoy 50 with the most seaward trajectory accelerated to SO cm s and 
Buoy 55, the inshore trajectory, accelerated to 145 envs. Ihese speeds are indicative of 
a strong narrow current along the northeast coast of Borneo with divergence to the 
south and north at the entrance to Macassar strait. This indicates a strong cyclonic 
shear was present across this current in that buoy separation in the Borneo Current was 
only about 50 km. 

From the AXBT survey conducted by the Naval Postgraduate School, the upper 
ocean temperature structure also shows the horizontal extent of these features. Hori- 
zontal depictions of the 23°C and the 14°C isotherms clearly identified the locations of 
the Mindanao and Halmahera eddies at the time of this survey. The most striking fea- 
ture observed from these data was the the frontal boundary between the SFC water and 
the NEC water as they retrollect then converge to form the origin of the NECC. Drifter 
speeds in excess of 120 cm, s were observed in this region at approximately 5°X, 
12S.5“E (Buoy 52). 

B. RECOMMENDATIONS 

Any conclusions draw’n from this study have limited scope and must be supple- 
mented with additional research. In this study there has been no correlation attempted 
between the buoy trajectories and the synoptic weather that occurred in June and July, 
19SS. This needs to be done. The seasonal evolution of the circulations associated with 
these eddies can only be inferred until this study is repeated during the northwest 
monsoon season. 

The Pacific-Indian Ocean tliroughflow is a stated objective of TOGA and WEPOCS. 
This study provided a glimpse of this throughllow via the trajectory of Buoy 55. How- 
ever, this is only a sample size of one. 1 highly recommend that any subsequent cruises 
to this region be extended to include a more thorough investigation of the mechanisms 
for the transport of north Pacific water to the Indian Ocean. 
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APPENDIX A. AXBT TEMPEKA I UKE PROI ILEb 



Data collection in support ol' the OXR sponsored Eddy Study occurred on the 
evenings of 11, 12, and 13 July 1988. Over the course of 33.2 air hours logged, 133 Air- 
borne Expendable Bathythermographs (AXBTs) were deployed, using channels 12, 14 
and 16. Reserve squadron VP-69 (from Whidby Island, WA) assets at NAS Cubi Point, 
RP provided the P-3 platform, as well as significant technical and logistical support. 
Additional coordination and support was provided by VP- 17 and the Naval 
Oceanography Command Facility at NAS Cubi Point. 

The intended collection procedure was to install an Airborne Digital Data Acquisi- 
tion System (ADDAS) on board the P-3 and digitally record and perform in situ analyses 
of AXBT measurements. The ADDAS was misrouted in shipping and did not arrive at 
NAS Cubi Point in time for scheduled use, precluding its use in the actual collection 
process. The data was recorded via on board BT profile output strip charts, anograms, 
and magnetic tape to be later transcribed on VllS tape using the ADDAS at the Cubi 
Point ASWOC, and post-processed at the Naval Postgraduate School. 

During collection, the average on station altitude was 2,500 feet and the average 
ground speed was 225 knots. There was approximately a 28% failure rate either due to 
AXBTs not "lighting olT', sending a weak signal, or due to a recording failure on board 
the aircraft. 

Included in this appendix are vertical profiles and listings of temperature at standard 
depths. The Longitudes gi\en are east. 
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TEMPERATURE (®C) TEMP 




(aa) 3dnss3bd 



53 



STATION: 4 LAT: 7 23.6 N LON: 128 9.7 
DATE: 7/12/88 TIME:1106Z 



TEMPERATURE (*C) PRESS TEMP 
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STATION; 6 LAT: 7 39.8 N LON: 128 50.9 
DATE: 7/12/88 TIME:1118Z 



TEMPCRATURE (®C) TEMP 
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STATION: 8 LAT: 7 14.9 N LON: 128 32.9 
DATE: 7/12/88 TIME: 1130Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION; 10 LAT: 7 14.4 N LON; 127 58.9 
DATE; 7/12/88 TIME:1141Z 



TEMPERATURE (*C) PRESS TEMP 
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STATION: 12 LAT: 7 15.1 N LON: 127 22.9 
DATE: 7/12/88 TIME:1153Z 



TEMPERATURE (*C) TEMP 
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STATION: 14 LAT: 6 49.8 N LON; 127 27.0 
DATE: 7/12/88 TIME:1211Z 



temperature: (*c) press 
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STATION: 16 LAT: 6 50.4 N LON: 128 1.0 
DATE: 7/12/88 TIME: 1236Z 



TEMPERATURE (®C) PRESS TEMP 
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STATIC*' LAT:6 50.2 N LON; 128 17.5 
DATE: 7/1^88 TIME: 1236Z 



irUPCRATURC (®C) PRESS TEMP 
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STATION: 18 LAT; 6 50.0 N LON: 128 34.0 V, 
DATE: 7/12/88 TIME:1241Z 



TEMPERATURE ("C) PRESS TEMP 
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STATION: 19 LAT: 6 50.0 N LON: 0 0.0 
DATE: 7/12/88 TIME: I248Z 



UMPFRATURF (“C) PRESS TEMP 





: o o */ 
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STATION: ?0 LAT; 6 ?/1.nN LON; 178 48.9 
fJATE: 7/12/88 TIME: 12532 



TEMPERATURE (®C) TEMP 
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STATION: 21 LAT: 6 25.0 N LON: 128 32.4 
DATE: 7/12/88 TIME:1300Z 



7FMPERATURE (*C) PRESS TEMP 
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STATION: 22 LAT: 6 25.1 N LON: 128 14.5 
DATE: 7/12/88 TIME:1306Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION; 23 LAT: 6 25.1 N LON; 127 56.2 
DATE: 7/12/88 TIME;1311Z 



TEMPERATURE (”C) TEMP 
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''TATION:24 LAT; 6 25.2 N LON: 127 33.5 -A' 
DATE: 7/12/88 TIME:1318Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION: 25 LAT: 6 25.1 N LON: 127 20.9 
DATE: 7/12/88 TIME: 1323Z 



TEMPERATURE (®C) 
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STATION: 26 LAT: 6 25.0 N LON: 0 0.0 
DATE: 7/12/88 TIME:1330Z 



TEMPERATURE (“C) TEMP 
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STATION: 28 LAT: 5 40.5 N LON: 125 54.1 
DATE: 7/12/88 TIME;1353Z 



TFMrCRATURE: (*C) PRESS TEMP 
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STATION: 29 LAT: S 43.4 N LON: 126 11.7 
DATE: 7/12/88 7IME;1400Z 



temperature: (®c) temp 
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STATION: 30 LAT: 5 48.1 N LON: 126 29.1 
DATE: 7/12/88 TIWE:1406Z 



TCMPFRATURE (*C) PRESS TEMP 
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STATION: 42 LAT: 5 46.8 N LON: 128 7.7 
DATE: 7/12/88 TIME: 1518Z 



TEMPCRATURt (®C) PRESS TEMP 
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STATION: 43 LAT: 5 41.7 N LON: 127 54.2 
DATE: 7/12/88 TIME:1530Z 



TCMPERATURE (“C) PRESS TEMP 
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STATION: 90 LAT: 6 13.9 N LON: 129 0.2 
DATE: 7/11/88 TIME:0853Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION: 91 LAT: 5 35.7 N LON: 129 0.1 
DATE: 7/11/88 TIME:0911Z 



temperature: (®c) press temp 
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STATION: 92 LAT: 5 21.0 N LON: 128 35.5 
DATE: 7/11/88 TIME:0918Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION: 93 LAT; 5 5.8 N LON: 128 34.9 
DATE: 7/11/88 TIME:0923Z 



7EMPFRATURE (“C) PRESS TEMP 
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STATION: 95 LAT; 4 12.0 N LON: 128 34.4 
DATE: 7/11/88 TIME:0936Z 



TEMPERATURE (*C) PRESS TEMP 
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STATION: 96 LAT: 3 33.0 N LON: 128 34.9 
DATE: 7/11/88 TIME:0948Z 



TEUPfRATURE (*C) PRESS TEMP 
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STATION: 98 LAT: 3 2.8 N LON: 128 35.8 
DATE; 7/11/88 TIME: 09532 



TEMPERATURE (*C) PRESS TEMP 
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STATION: 99 LAT: 2 23.8 N LON: 129 0.2 
DATE: 7/11/88 T1ME:1000Z 



TEMPERATURr (*C) TEMP 
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STATION: 100 LAT: 2 8.5 N LON: 129 4.8 
DATE: 7/11/88 TIME:1011Z 



TEMPERATURE (“C) PRESS TEMP 
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STATION: 102 LAT: 1 14.8 N LON; 129 2.5 
DATE: 7/11/88 TIME:1023Z 



TCk^PrRATURE: (®C) PRESS TEMP 
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STATION: 108 LAT: 3 20.0 N LON: 129 14.5 
DATE: 7/11/88 TIME:1053Z 



TCMPERATURE ("C) PRESS TEMP 
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STATION: 109 LAT: 3 32.8 N LON: 129 14.6 
DATE: 7/11/88 TIME:1106Z 



TEMPERATURf (®C) PRESS TEMP 
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STATION; 110 LAT: 4 12.1 N LON: 129 13.5 
DATE; 7/11/88 TIME:1111Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION: 111 LAT: 4 27.1 N LON; 129 15.7 
DATE; 7/11/88 TIME: 1118Z 



TEMPfRATURE (*C) PRESS TEMP 
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STATION: 112 LAT: 5 5.8 N LON: 129 15.8 
DATE: 7/11/88 TIME: 1123Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION- 113 LAT:5 21.0 N LON: 129 16.3 
DATE:;, i 1/88 TIME:1130Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION: 114 LAT: 6 0.1 N LON: 129 15.8 
DATE: 7/11/88 TIME; 1136Z 



TEMPERATURE (*C) PRESS TEMP 
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STATION: 117 LAT; 6 1.6 N LON: 129 29.5 
DATE: 7/11/88 TIME:1148Z 



TEMPERATURE (“C) TEMP 
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STATION; 118 LAT: 5 21.0 N LON; 129 29.1 
DATE: 7/11/88 TIME:1200Z 



TEMPERATURE (“C) PRESS TEMP 
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STATION: 119 LAT: 5 5.8 N LON: 129 27.1 
DATE: 7/11/88 TIME:I206Z 



TEMPERATURE (®C) TEMP 




(ea) 3unss3ad 



95 



STATION: 120 LAT: 4 27.0 N LON: 129 29.0 
DATE: 7/11/88 TIME:1211Z 



TEMPERATURE (®C) PRESS TEMP 
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STATION: 122 LAT: 3 33.0 N LON: 129 29.5 
DATE: 97/11/88 TIME:1223Z 



TEMPFRATURE (®C) PRESS TEMP 



o lO o 
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STATION: 123 LAT: 3 18.0 N LON: 129 25.4 
DATE: 7/11/88 TIME: 1230Z 



TEMPERATURE (®C) TEMP 
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STATION: 124 LAT: 3 2.8 N LON: 129 30.0 
DATE: 7/11/88 T1ME:1241Z 



TrWPERATURE (“C) PR^SS TEMP 
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STATION: 126 LAT; 2 27.4 N LON: 0 0.0 
DATE: 7/11/88 T1ME:1248Z 



TEMPERATURE (*C) PRESS TEMP 
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STATION: 127 LAT: 1 30.4 N LON: 129 30.0 
DATE: 7/11/88 TIME: 1253Z 



temperature: (®c) press temp 
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STATION: 128 LAT: 1 15.3 N LON: 130 9.1 
DATE: 7/11/88 TIME:1306Z 



APPENDIX B. LAGIUNGIAN TEMPERATURE I IME SERIES 

This appendix is included to provide the temperature time series which were dis- 
cussed in this study and may be of interest to the reader. Also included are the temper- 
ature time series based on the raw data. 
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Sea Surface Temperature, Buoy 51 
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Sea Surface Temperature, Buoy 52 
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Sea Surface Temperature, Buoy 53 
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Sea Surface TemporaturQ, Buoy 54 
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SurfacQ Tomportature, Buoy 55 




s 




108 



28.5 




LIST OF REFERENCES 



Emer}', W.J., W.G. Lee, and L. Magaard, Geographic and Seasonal Distributions of 
Bruni-Vaisala frequency and Rossby Radii in the North Pacific and North Atlantic 
J. Phys. Ocean, 14, 294-317, 1984 

Cannon, Glenn A., Characteristics of the Waters East of Mindanao, Philippines Islands, 
August 1965, in The Kuroshio, pp. 205-211. 

Lukas, R., Firing, E., Moore, D., Helsley, C. Research Proposal Submitted to the 
Physical Oceanography Section of the National Science ¥ oundaiion, Low- Laiiiude 
IVestern Boundary Currents in the Pacific Ocean, 1987. 

Niiler, P.P., Davis, R.E., White, M.J., Water-following Characteristics of a Mixed Layer 
Drifter, Deep-Sea Res, 34, 1S67-1881, 1987. 

Poulain, P.M., Illeman, J.D., Niiler, P.P., Drifter Observations in the California Current 
System (1985,1986), University of California Scripts Institution of Oceanography 
Data Report. 1987. 

Takahashi, T., Hydrographical Researches in the Western Equatorial Pacific, Memoirs 
of the Faculty of Fisheries Kagoshima University 7, 141-147, 1959. 

Tsuchiya, M., Lukas, R., Western Equatorial Pacific Ocean Circulation Stud y 
(WEPOCS 11) Shipboard Chemical and Physical Data Report, NSF, 1987 

Wrtki, K., Magaard, L. and Hager, J., Eddy Energy in the Oceans J. Geophyi. Res., 81, 
2641-2646, 1976. 

Wyrtki, K. AAGA Report, Vol. 2: Scientific Results of Marine Investigations of the 
South China Sea and the Gulf of Thailand, 1961 



109 



INITIAL DISTRIBUTION LIST 


No. Copies 


1. Defense Technical Information Center 
Cameron Station 
Alexandria, VA 22304-6145 


2 


2. Librar>’, Code 0142 

Naval Postgraduate School 
Monterey, CA 93943-5002 

1. Chairman (Code 63Rd) 

Department of Meteorology 
Naval Postgraduate School 
Monterey. CA 93943-5000 

2. Chairman (Code 68Co) 

Department of Oceanography 
Naval Postgraduate School 
Monterey, CA 93943-5000 

3. Professor C. A. Collins(Code 68 Co) 
Department of Oceanography 
Naval Postgraduate School 
Monterey, CA 93943-5OO0 

4. Professor Mar>- Batteen (Code 68BV) 
Department of Oceanography 
Naval Postgraduate School 
Monterey, CA 93943-5000 

5. Director Naval Oceanography Division 
Naval Observatory 

34th and Massachusetts Avenue NW 
Washington, DC 20390 

6. Commander Naval Oceanography Command 
Naval Oceanography Command 

Stennis Space Center, MS 
39529-5000 

7. Commanding Officer 

Fleet Numerical Oceanography Center 
Monterey, CA 93943 

8. Commanding Ofl'icer 
Naval Oceanographic Office 
Stennis Space Center, MS 39522-5001 


2 

1 

1 

1 

1 

1 

1 

1 

1 



110 



9. Commanding Officer 

Naval Ocean Research and Development Activity 
Stennis Space Center, MS 39522-5001 

10. Commanding Olficer 

Naval Environmental Prediction Research Facility 
Monterey. CA 93943-5006 

11. Chairman, Oceanography Department 
U. S. Naval Academy 

Annapolis, MD 21402 

12. Chief of Naval Research 
800 North Quincy Street 
ArUngton, VA 22217 

13. Office of Naval Research (Code 420) 

Naval Ocean Research and Development Activity 
800 North Quincy Street 
Arlington, VA 22217 

14. LCDR Glen H. Carpenter 
4266 Greensboro Rd 
Corpus Christi, TX 7S413 

15. Scientific Liaison Office 
Ollice of Naval Research 
Scripps Imitution of Oceanographv 
La Jolla, CA 92037 

16. Library 

Department of Oceanography 
University of Washington 
Seattle, \VA 98105 

17. Dr. John Toole 

Physical Oceanography Department 
Woods Hole OceanoEraphic Institute 
Woods Hole .MA. 02543 

18. Research /idmiiiistration 
Code 012 

Naval Postgraduate School 
Monterey, CA 93943 



19. Professor Klaus Wyrtki 
Department of Oceanography 
Universitv of Hawaii, Manoa 
Honolulu. HI 96822 

20. Dr.Eric Fleming 
Department of Oceanography 
University of Hawaii, .Manoa 
Honolulu, HI 96822 



21. Dr.Roger Lukas 
Department of Oceanography 
University of Hawaii, Manoa 
Honolulu, HI 96822 

22. Dr. Peter Hacker 
Department of Oceanography 
University of Hawaii, Manoa 
Honolulu, HI 96822 

23. Dr. Philip Richardson 

Physical Oceanography Department 
Woods Hole Oceanographic Institute 
Woods Hole MA. 02543 

24. Dr. Rana Fine (Code 68BV) 

RSxMAS 

University of Miami 

4900 Rockenbacker Causeway 

Miami, FL 33149 



■pTTT^ijW.V TT'TOX LlBBi-’' 



Thesis 
' C2635 
c.l 



Carpenter 



Thesis 

C2635 Carpenter 
c.l Surface circulation 

associated with the 
Mindanao and Halmahera 
Eddies . 



thesC2635 

Surface circulation associated with the 




3 2768 000 82690 3 

DUDLEY KNOX LIBRARY 



